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We model a small quantum dot with a magnetic impurity by the Anderson Hamiltonian with a 
supplementary exchange interaction term. The transport calculations are performed by means of 
the Green functions within the equation of motion scheme, in which two decoupling procedures are 
proposed, for high and low temperatures, respectively. The paper focuses on the charge fluctuations 
for such a system, aspect not addressed before, as well as on the Kondo resonance. We show a 
specific role of the excited state, which can be observed in transport and in spin-spin correlations. 
^ ' Our studies show on a new many-body feature of the phase shift of transmitted electrons, which is 

manifested in a specific dip. In the Kondo regime, our calculations complement existing theoretical 
results. The system shows three Kondo peaks in the density of states: one at the Fermi energy and 
two side peaks, at a distance corresponding to the singlet-triplet level spacing. The existence of the 
central peak is conditioned by a degenerate state (the triplet) below the Fermi energy. 
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I. INTRODUCTION 



. Quantum dots (QD) arc promising systems for nano-technology because their quantum-mechanical parameters are 
O ' easy to control. Motivated, in part, by the recent developing in the field of magnetic semiconductors or by the 
, O, . spintronic advances, theorists and experimentalists have started to address also the problem of quantum dots with 
magnetic impurities @j S S S IB- l3] ■ In such systems, electrons in the quantum dot are subject to a supplementary 
P\| i exchange interaction. This induces a competition between the formation of the Kondo cloud and the magnetic 
J> • interactions inside QD, but also specific charge fluctuations, aspect less addressed in literature. Sometimes, the role 
\ of the impurity is played by another quantum dot with odd occupancy @, but also, as in a recent experiment the 
1^ I impurity can be a transitional magnetic atom (Co). The splitting of the Kondo peak was attributed by the authors 
I to the indirect RKKY interaction & [l^, [ll| . We believe that a central peak would have also been present for a 
, much lower temperature. Govorov [2l studied theoretically the spectrum and the response to optical excitations of an 
\^ ■ InGaAs / GaAs quantum dot with a single Mn impurity. Recently, Murthy Q also addressed the problem of interplay 
, between the Stoner and the Kondo effect in a large QD with a magnetic impurity. 

■ The competition between the Kondo and antiferromagnetic couplin g w as also intensively studied in a double-QD 
^ ; geometry coupled by hopping or by the RKKY interaction (e.g. [1, [fills', 'T^, Tsl, [13 ) . This is related 
H ■ with another interesting issue: the Kondo effect in QD with an integer spin [1^, i20, 21, 23, [2^. In the experiment of 
] Sasaki [2^, the applied magnetic field can change the ground state and the relative position of the singlet and triplet 
states, and it can also induce a degeneracy point in the energy spectrum. In the vicinity of the degeneracy point, 
three Kondo peaks are noticed in the differential conductance, which is in agreement with our calculations for the 
^ Kondo regime. 

;V ■ Previous theoretical papers have focused mainly on the Kondo effect in large QDs with magnetic impurities, or on 
the two-impurity problem. In the present paper we are also interested in the Kondo resonance, but in a different 
system: a small QD with a magnetic impurity that is coupled only to the electrons in the dot. We are interested in 
5— i ' interplay between magnetic interactions: the dot-leads Kondo coupling (as in the experimental set-up [2^. [25j) and 
the intra-dot magnetic exchange interactions (between electrons in QD and the magnetic impurity). We would like 
to consider spin-flip processes on the impurity and their detection by the current flowing through the QD, both in 
high and low temperature regimes. In spectral terms, the flip of the impurity spin is associated with singlet-triplet 
excitations. To our knowledge, the problem was not considered yet. 

Apart from the Kondo regime, we shall also analyze the regime of high temperatures. The solution in this case is 
more easy to obtain, and it allows simple analytical expressions. These will be used to discuss the charge fluctuations 
and spin-spin correlations aspects. An electron that enters QD forms with the magnetic impurity a correlated singlet 
or a triplet ground state, respectively to the sign of the exchange parameter. Both situations are met in experiments 
0. In these cases, the excited states also play an important role. It can be seen in the conductance peaks, which 



FIG. 1: Scheme of a quantum dot coupled to ideal leads through tunnelling barriers t. A magnetic impurity interacts with 
electrons in the dot through an exchange coupling J. 

amplitude strongly depend on the coupling to the leads, and in the spin-spin correlation functions, especially when 
the excited states lie below the Fermi level. We also address the problem of the transmission pha se. The measurement 
of the transmission phase (and not only the amplitude) is possible to realize, for some time j26i [27|. and can provide 
interesting supplementary information. For our system, a dip can be noticed in the phase evolution, precisely at the 
excited state position, in condition of constant electronic occupancy on QD. We will discuss our results in connection 
with a recent experiment [28j . where similar features were observed. Our studies are an explicit demonstration of a 
generalized Fridel sum rule in the presence of electronic correlations, the theorem which was just recently proved by 
Rontani [H. 

The paper is organized as follows: in the Section II and III we describe our model and the Green functions approach. 
The equation of motion (EOM) technique is presented in detail. A decoupling scheme adequate for high temperatures 
is presented in section IV. The scheme allows to obtain analytical results for the conductance, dot occupancy, spin-spin 
correlator and the phase evolution. In Section V, we get an insight into the Kondo regime using a more advanced 
decoupling scheme. One can observe the formation of three Kondo peaks in the local DOS: at the Fermi energy and 
two side peaks. The central peak exists only if a degenerate state (the triplet) is below the Fermi energy. Section VI 
contains concluding remarks. For completeness, the spectrum of the system and a discussion on the formation of the 
Kondo resonance are presented in the Appendix A and B, respectively. 

II. MODEL OF A QUANTUM DOT WITH A MAGNETIC IMPURITY 

The considered system of the QD with a magnetic impurity is presented in Fig. 1. Conduction electrons passing 
through the QD interact with accumulated electrons leading to the Kondo resonance in low temperatures. A magnetic 
impurity is connected to the electrons in the QD by the exchange coupling J. The Hamiltonian corresponding to our 
model is written as 

-ff= ^kcl^^aCka^a + f^Q^cl^Coa + Ucl^Co'^cl^Coi + Js- S + ^ ta{cl„Cka,a + h.C.) , (1) 

where the first term represents the electrons in the lead a = L, R; the second term stands for the up-most electronic 
level eo in the dot. The the third and the fourth terms describe the interactions: the Coulomb interaction of electrons 
with opposite spin orientation at the level eo and exchange interactions with the magnetic impurity. Here, s is the 
spin operator for electrons at the level eg, whereas S denotes the spin operator at the magnetic impurity. It is imposed 
a single occupancy at the impurity -t- d^^di = 1 and the the spin S = 1/2. Therefore, = l/2(d|(if — djdx), 

5+ = d''^ di, = rfjrfti ^^'^ the exchange term between the impurity and the spin of the electron localized at the QD 
can be written as 

Js - S ^ ^^(cJtCot - cJj^cox)5"^ + ^Jcl^coiS~ + ijcJj^co|S'+ . (2) 

The last term in the Hamiltonian (1) corresponds to the coupling between the quantum dot and the leads. We are 
interested in specific correlations effects, so the eventual effects introduced by an asymmetric coupling are outside our 
purpose here. It shall be considered, for simplicity, ti, = ta = t. 

The energy spectrum of the isolated QD is easily to calculate (see the Appendix A). We will show that the ground 
states, as well as the excited states, participate in the electronic transport. The conductance peaks corresponding to 
the transport through the excites states are smaller, but they are high enough to be observable in an experiment. 
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Our model (1) is a generalization of the single impurity Anderson model (SIAM). to which a spin- interaction term 
is added. Recently, Zhang et al. [l^l considered also an extended SIAM model with a a one-body phenomenological 
spin-flip term. Their purpose was to study a different situation, namely, the electronic transport through a quantum 
dot with spin-orbit interactions coupled to ferromagnetic leads. In the Hamiltonian (1), the total spin is conserved 
(in contrast to the model used in [30|), which is important for the electronic transport. 



III. DETERMINATION OF GREEN FUNCTIONS IN EQUATION OF MOTION SCHEME 

Now, we want to determine the conductance and other physical quantities for the model described by the Hamil- 
tonian (1). The current can be expressed by means of the non-equilibrium Green functions as [Sli, [32| 

J^^^l dcjihicj) - /fl(^)](-Im((coTl4))), (3) 

where F = 2nt^p and p ~ 1/2D is the DOS for the square band approximation; the half-width D will be taken 
as unity. Jl and fa are the Fermi distribution functions in the left and right leads, respectively, ((cof |cj|)) is the 
retarded single particle Green function for an electron with the spin a =1 at the QD, which can be determined by 
the equation of motion (EOM). The differential conductance for the quasi-equilibrium case can be deduced from the 
formula (3) 

S-^-^rJ M-^)lm{{co,\cl^)) . (4) 

The equation of motion (EOM) for the energy dependent retarded Green function is given by 

c.{{Am = {{A,B}) + {{[A,H]\B)) , (5) 

where {{A, B}) is the thermal average of the anticommutator between the operators A and B. Whenever necessary, 
the averages will be calculated by means of the fluctuation-dissipation theorem {AB) = — (I/tt) J dujf{u!)lm{{B\A)) , 
where /{tij) is the distribution function at the dot. In equilibrium, it is equal with the Fermi distribution in the leads. 
To determine the current (3) one needs the single-particle Green fimction ((cocr|cJ^)), for which the EOM is 

{a; - eo){{co^\cl^)) = 1 + U{{co^cl^coi\cl^)) + ^J{{co^S'\cl^)) (6) 

+ ^^((coi5-|cS^))+2t^((cfeTl4» 

fc 

As usual for an interacting problem, one has also to calculate the many-particle Green functions 

(a;-6o)((coT5^|ct^)) = (5^> + i J((cot kj^)) + [/((coTcJ^coi^lct^)) (7) 

k 

{u;~eo + ^J){{coiS~\cl^)) = ^ J((co||cJ^)) + ([/ + ^J){{cqic1^co^S-\cI^)) 

- Ij{{co^S^\cI^)) + J{{co^clcoiS^\cl^)) 
+ 2tJ2{{ckiS-\cl^)) 

k 

(cj - eo - C/)((cotcJ^cox|cJ^)) = (cJ^cqx) + ^J{{co^cl^CoiS'\cl^)) + ^ J((coicJ^CoT5'" |cj^)) 

+ 2t^((cfc|cJ^coi|C|!,^)) - 2t^((coT4^coi|cJ^)) 
fc fc 

+ 2t^((cotcJ^Cfcx|C|!,^)) 
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and 

(w-eo-C/)((coTcJ^coiS'^|cJ-p)) = (cJ^cqxS'^) + ^ J((cotcJ^coxl4^)) + ^ J((coi4^coTS'~|cJ-p)) (8) 

k k 

+ ^tY,{{coAiCkiS'\cl^)) 

k 

(w- eo - - U){{coic1^cotS-\cI^)) = -(cJ^cox5") + ^J{{co^cl^coi\cl^-)) + ^J{{co^cl^coiS''\cl^-)) 

+ 2t^((cfcicJ^cotS'"|cJ^)) - 2t^((co44^coT5'"|cJ^)) 

fc 

k 

In the absence of a magnetic field, one can consider (S^) = 0, (?^oi) = ('^ot) = (co|Cot) and (cJj^co^S'^) = 
— i(cJ|Cox5^) = ^■|(** ■ 'S'): ("S*^)^ = J (see also [S^ for some details on how the constriction (i|d| + djdi = 1 is 
used) . 

We have got a set of six equations, Eqs.(6)-(8), for all Green functions at the quantum dot. These functions are 
coupled with other many-particle Green functions {{cka A\cq^)) , which describes the coupling between QD and the 
electrodes. The EOM procedure requires to write also the equation of motion for the new set of Green functions, which 
will be coupled with another new set of the higher order many-particle Green functions. In order to close the infinite 
series of the equations of motion, one has to decouple many-particle Green functions. The simplest approximation is 

2tY,{{c,^,A\clj) « -^^{{co..A\cl,)) . (9) 

k 

The approximation takes into account charge fluctuations between the electrode and QD, but it neglects the spin- flip 
processes and the Kondo resonance. This approximation correctly describes the limit of the Coulomb blockade and 
gives the same results as those obtained by Nagaoka [s^ and by Meir et al [s^l (for QD without impurity) in high 
temperatures. It is important to notice that all Green functions at QD are treated exactly, so the intra-dot correlations 
are fully accounted. This approximation will be used in the next section to study the high temperature properties of 
our model. 



IV. STUDIES OF THE HIGH TEMPERATURES REGIME 

Using the approximation ([9]) we get a set of six linear equations, with the correlators (?io|) and (cJ|Co|S'^) as 
parameters. One can easily find analytical expressions for all many-particle Green functions at the QD. Next, the 
correlators arc calculated by means of the appropriate Green functions (no|) = — (I/tt) J (iti>/(w)Im((co| |cJ|)) and 
(cJ|-Co|5'^) = — (I/tt) j rfaj/(w)Im((co|S'^|cJ-|.)). In order to calculate these integrals one has first to decompose the 
Green functions in simple fractions. The solutions are expressed by the digamma function [ssj . which at zero tem- 
perature, turns into the arctan function. 

The main results of this section are presented in Fig. 2. The upper panels show the conductance (solid curves) 
and the electron occupancy of the QD (dotted curves) as a function of the position of the dot energy eo for different 
parameters T. The lower panels show the corresponding behavior of the spin-spin correlator. 

For J > - column a) - four conductance peaks are clearly seen, while for J < - column b) - there are two 
pronounced peaks and two much smaller peaks. To explain the origin of the peaks, let us start with a high value 
of eoj which is much above the Fermi energy Ep ~ 0. For this case the dot is empty. When eg is lowered (by an 
applied gate potential), a first electron is allowed to jump into the QD on the lowest available level. This is the singlet 
state, for case a) -see Appendix A- and it is located at cq = 0.15. For case b), the lowest state is the triplet state at 
eo = 0.05. By further lowering of cq, the excited state energy becomes aligned with the Fermi level. It is the triplet 
state at eo = —0.05, or the singlet state at eo = —0.15, for the case a) and b) respectively. At these positions the 
conductance shows peaks, for which the height is smaller than for those corresponding to the ground states. We do 
not observe any changes of the electron occupancy at these positions. 

The transmission through the excited states, where the total occupancy on the dot is constant, can be explained 
by the overlapping with the ground state and fluctuations of the charge distribution between levels. The electron 
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has a probability to exist on the excited state, and from here it can tunnel into the leads. This explains the strong 
coupling-dependence of the excited states conductance. Notice a significant increase of the peaks corresponding to 
the transport through excited states as the coupling F increases from 0.01 to 0.08. Experimentalists indeed have 
measured such small conductance peaks corresponding to transport through excited states (e.g. (36l. \3^). when the 
total occupancy on the QD is constant. Apart from describing an impurity dot, our model may be one of the simplest 
to allow the analysis of excited states behavior in transport. 

The transmittancc through the excited states is more pronounced in Fig2.a than in Fig. 2b. This is because the triplet 
state (that is the excited state in the case a) is three-fold degenerate and in favored in transport in comparison with 
the non-degenerate singlet state. From the same reason, one can notice that, on the contrary, the peak corresponding 
to the ground state (triplet) in Fig. 2b is higher than in Fig. 2a (singlet). The conductance peaks corresponding to 
the transport through the ground states remain then practically unchanged, only a broadening effect is noticed. A 
detailed discussion of the peaks amplitude will be given further in this section. 

We remind that the resonances in conductance appear when two quantum states with consecutive occupancy are 
energetically available in the same time, allowing for charge fluctuations. So, one can find the positions of the 
resonances with the equation E(iy,n) = E{v' ,n ± 1), meaning that an energy from the spectrum with n electrons 
becomes equal with an energy for n ± 1 electrons - see Appendix A. By further lowering of eg, the energy for double 
occupancy will become equal with the excited single-particle energy. This will generate a small peak in conductance, 
and no variation of occupancy. Finally, a big peak is noticed corresponding to the adding of the second electron. 
Actually, a particle- whole symmetry is easy to notice in Fig. 2, and we shall focus the further discussion on the first 
two peaks from the right (marked by the vertical, dotted lines). 

The lower panel in Fig. 2 shows the evolution of the spin-spin correlator (s • S), which is a measure of the coupling 
between the spin of an electron accumulated at the QD and the magnetic impurity spin. For J > 0, electrons at 
the QD are antiferromagnetic coupled with the impurity spin, and therefore, the spin-spin correlator is negative. Its 
minimal value is —0.7, —0.5 and —0.3 for F = 0.01, 0.04 and 0.08, respectively. For the case F = 0.01, the length of 
the total spin, {S^^^ — (s^) + 2(s • S) + (5^) is then reduced to 0.1, i.e. it is almost compensated. When the triplet 
state gets below Fermi energy, the correlator (s - S) increases, because the triplet state favors ferromagnetic coupling. 
Its maximal value is < s • 5* >« —0.1, which means < 3"^^^ >w 1.3 close to 3/2 - the value of the magnetic moment for 
two free 5* = 1/2 spins - see the central part of the graph. The result is important because it shows that the excited 
states contribute to the formation of the local magnetic moment at the QD. It can have also influence on the Kondo 
singlet formation with the electrons from the leads. Section V will show that indeed this is the case. 

The evolution of the spin-spin correlator for the ferromagnetic case, presented in Fig. 2b, can be explained in the 
same way. When the triplet state crosses Ep the spin-spin correlator reaches maximal values 0.23, 0.2 and 0.15, for 
F = 0.01, 0.04 and 0.08, respectively. The total magnetic moment increases up to 1.94, 1.86 and 1.7. For each of the 
three situations, the spin-spin correlator and the total moment begin to decrease when the singlet state goes below 
Ep . 

We want to discuss in more details the heights of the conductance peaks corresponding to the singlet and the triplet 
states. A simple analytical solution can be obtained for the conductance, for T = and in the large U limit , for the 
singlet and triplet resonances 

c-c +r 1 + p 3(l + 2</.s) 

Q-Ss + Gt- 2i3 + ^s-^scl>Tf' + 2(3 + 40s-40s0t)''^ ' ^''^ 

where Rs = (2eV/i) x rVp^ + {Ep - eo + 3J/4)2] and Rt = {2e^/h) x V^/[T^ + [Ep - eo - J/'^f] corresponds to the 
resonant conductance through the singlet and the triplet level, respectively, (ps = arctan[(eo — Ep — 3J/4)/r]/7r and 
(pT = arctan[(eo — i?_F + J/4)/F]/7r. The fractions that multiply the resonances Rs{t) contain the information about the 
height of the conductance peaks. The ground and excited states role changes with the sign of J, and the dependence 
on the coupling also changes, being much more pronounced for the excited state. All these effects are incorporated in 
formula (10). This zero-temperature formula is given for its simplicity, but it describes well also situations like in Fig. 2, 
where temperature T = 10~^ is considerably lower than the broadening due to coupling with the leads (F = 0.01 to 
0.08). For higher temperatures, one has to replace in formula (10): 4>s = — Im{vI'[l/2— i(eo— i5_F— 3J/4-|-ir)/(27rr)]}/7r, 
(j)T in a similar way is the digamma function, see also (38jl. and then perform an integral required by the conductance 
formula deduced from Eq.3. Fig. 2 was plotted in this way, with the complete formula. 

By employing further the simplified form Eq.(lO), one can estimate the ratio between the heights of the conductance 
for the singlet and triplet peaks. We notice that in Eq.(lO), the parameter (ps (0r) vanishes at the singlet (triplet) 
resonant level, and therefore, one can straightforward calculate the ratio of the conductance peaks 

max[gs] ^ (l + 20o)(3-40o) 

maxie^T] 9(1-200) ' ^ ' 
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FIG. 2: The conductance (solid curves), the dot occupancy per spin (noo-) (dotted curves) - upper panels, and the spin-spin 
correlator {s-S) - lower panels, plotted versus the dot energy eo, for U = 0.5, |J| — 0.2, Ef — 0, T = 10"'^, and for different 
values of r = 0.01, 0.04 and 0.08. The vertical lines indicate the position of the Singlet and the Triplet levels. 



where 4>o = arctan( J/r)/7r. For the antiferromagnetic coupling the electronic transmission through the singlet state 
dominates and in the hmit J ^ F this ratio goes to infinity and the triplet peak disappears. In the case of the 
ferromagnetic coupling and | J| 3> F we have opposite situation - the singlet states disappears from transport and 
only the peak corresponding to the triplet ground state is visible. 

In 0, [3^, the authors model a similar spin-dependent transport situation (for the QD and QPC, respectively), but 
the correlations and the effects of the Fermi sea are neglected. Their problem is analogous to the electron collision 
with a hydrogen atom (see In such a case, by averaging all possible scattering processes, the height of the 

singlet and triplet peaks will be proportional with their degeneracy resulting the 1/3 peak height ratio - independent 
of the sign of J, also independent of the coupling to the leads. Formulas (lO)-(ll) are a natural generalization. The 
difference results from the charge accumulation and the Coulomb blockade effect. 

In the final part of this section, we analyze the evolution of the phase shift for the scattered electrons. The phase 
problem has raised considerable interest in mesoscopic physics [4l[. The principle of the phase shift experimental 
measurement p6l [27j uses an interference process, which has been realized by inserting QD - for which the phase is 
measured - in one arm of an Aharonov-Bohm interferometer. The electronic wave passing through QD acquires a 
phase shift and then interferes with the wave travelling through the reference arm. In this way, one can measure not 
only the amplitude of the transmission, but also its phase, which is extracted from the phase of the Aharonov-Bohm 
conductance oscillations. The authors expected a phase evolution in concordance with the Fridel sum rule, growing 
with TT on resonances, and with constant values between them. Instead, an intriguing universal phase-lapse effect 
was noticed [2^ H^, meaning that the phase has a drop of tt between any pair of resonances. The aspect has not 
yet received a satisfactory explanation. Recently, Avinun-Kalish et al. |28| realized a similar experiment, but with 
a supplementary possibility to control the exact number of electrons in QD (up to 20 electrons). For the first few 
electrons (about 10), the phase shift had a complex non-universal behavior, with dips rather than phase-lapses, and 
for the next electrons, the "universal" behavior was recovered. Rontani demonstrated, in a very recent paper [2^, 
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that such features in the phase evolution are normal to happen, as the mesoscopic systems with electronic correlations 
do not obey the classical Fridel sum rule. It means that the experimentally measured phase shift is not equal to tt 
multiplied with the accumulated charge, but may have variations in the case of electronic correlations in the system. 
We explicitly prove the existence of such dips in the phase evolution, for the model discussed in this paper. Moreover, 
they receive a spectral interpretation, as being a fingerprint of the excited states. 

For simplicity, we consider here only the situation F = 0.01, when the singlet and triplet peaks are well-separated 
(but the conclusions are general). One could naively expect that the phase grows with 2tt on the two separate peaks, 
or that the phase evolution obeys the Fridel sum rule and accurately follows the occupancy curves. However, both 
hypotheses are wrong, as the ground state is decisively influencing the phase evolution on the excited state. We must 
remember that the very existence of the excited states peaks is attributed to the overlapping with the ground states 
peaks. 

The phase shift can be expressed as the argument of the single-electron Green function 0, H^l 



$ = arctan 



Im((co^|4^))^^j 
Re((coo-|4^))i^=i 



(12) 



In a single-particle scattering picture, the phase defined as (12) is equal to the phase-shift of the electron wave- 
function. The definition is identical with the definition of 9t in [i^]. The authors define the phases as the argument 
of the scattering matrix elements, which are the retarded Green functions. 



71/2 




FIG. 3: Evolution of the phase shift $ (solid curve), the real and the imaginary part of the single particle Green function 
G = {{coa-\cli^))^=Ep (dashed and dotted curves - in 1/D units) for the same parameters as in Fig. 2a and F = 0.01. 



Fig. 3 presents the phase evolution, the real and the imaginary part of the Green function. The phase evolves with tt 
on the peak corresponding to the ground state (around eq = 0.15), where an electron is added. Then the phase presents 
a dip on the excited state peak (around eq = —0.05), in conditions of unchanged dot occupancy (see Fig. 2). Notice 
that the real part of the Green function (the dashed curve in Fig. 3) changes sign at the ground state position, but 
not at the excited state position, which is why the argument of the complex Green function (12) remains in the same 
quadrant in the latter case. The phase exactly at the triplet position is easy to calculate. We obtained the minimal 
value in the dip (tt — arctan[6/7r]) ~ O.GStt |43| . for a low coupling (see Fig. 3). For the negative exchange parameter 
(as in Fig. 2, column b), and the singlet excited state, the minimal value for the dip is (tt — arctan[2/(37r)]) ~ 0.937r. 

The observed dip is a signature of the excited states and can be used to probe the existence of such states through 
phase measurements. Also, phase measurements can, in principle, provide information about the nature of the excited 
states - we obtained the amplitude of the dip different for the singlet and triplet. 

Our phase calculations can be correlated with the recent experiment by Avinun et al. p8| . where such dips in the 
phase evolution were observed. As mentioned, the authors measured the transmission phase, and have reported a 
non-universal behavior. The measured phase showed the expected evolution with tt on the resonances, but the authors 
also observed a number of dips of variable amplitude (but less than tt). Rontani (29j proved that the phase shift is 
not necessarily an integer of tt between two consecutive Coulomb blockade regions, if a complete formula is employed. 
The author stated that values lower than tt for the phase, through regions with constant occupancy, are a fingerprint 
of electronic correlations. The dips may even reach the value zero for the case of the spin blockade. 



V. AN INSIGHT ON THE KONDO REGIME 



In this section, we analyze the Kondo physics for our system, in a low temperature. In this regime, the electron 
from the dot may develop a Kondo-corrclated state with the electrons in the leads. This is in competition with the 
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spin-coupling with the magnetic impurity. 

A study of the Kondo resonance requires a method which takes into account correlations with electrons in the leads 
and the virtual spin-flip processes. One has to use the higher order decoupling scheme than that one used in the 
previous section. The next step of EOM involves the Green functions with one operator from the leads and decoupling 
for the Green functions with two operators from the leads. Here, an example of EOM 

q 

+ \j{{ckicl^co'^S-\cl^)) - ^J{{ckicl^coicl^coiS-\cl^)) . (13) 
The higher order Green functions are decoupled according the following scheme 

{{co]cl^Cki\cl^)) w (4iCfci)((coT|cJ|)) ~ 4g/(efc)((coT|c5|)) , (14) 

{{co^clckiS^lcl^)) « 4J(efe)((coT^1ct^» + 4,A/(e,)(cVoi^')((coT|cJT» • 

The decouphng procedure follows that one proposed by Meirjs^l (also used e.g. in [10, H, HE]) and by Nagaoka (s^l 
for the Green functions that also contain the localized spin operator. Only the averages that conserve the total spin 
are kept. The system of equations must be closed, so we need to express further the spin-spin correlator (cJ|CfcxS'~) 
between the electrons in the leads and the localized spin. This is considered to be proportional to the dot spin-spin 
correlator, by a parameter A < 1. The choice is justified by the fact that the localized spin is only indirectly coupled 
to the electrons in the leads. We have checked that the main results are not qualitatively influence for any < A < 1. 
In the numerical calculations we put A = t'^/D^. 

After the decoupling scheme (14), the set of Eqs.(13) (there are 14 equations) becomes linear and is easy to solve, 
in the U oo limit. The solution for the Green functions contains the fractions: 1/{lo — e^), 1/{lu — ej, + J) and 
1/{lu — efc — J). Next, we come back to the initial system Eqs.(6)-(8), where the summation over k leads to the 
appearance of three digamma functions [sHi 

^ / ^ /(ffc) TT , /27rT\ ^ /I uj + a\ 

with = and ±J. Here, ^ denotes the digamma function. For T < Tk, the function Fa has a pronounced maximum 
around lo = —a. 

Fig. 4 presents the density of states (DOS) of our system for three different cases: a) when both levels, singlet 
and triplet, are below Ep, b) the triplet state is below Ep, but the singlet state is above, and c) the singlet state 
is below Ep and the triplet above Ep. In the panel a), the singlet and triplet peaks overlap, due to the small J /T 
ratio. Our model gives three peaks in the DOS: one at the Fermi level Ep = and two symmetrical side peaks, at 
±J, corresponding to the basic excitations in the system. The central peak can be associated with the three-fold 
degeneracy of the triplet state. We point out that in the measurements of Sasaki (see Fig. 4a in [11]) such a three-peak 
structure can be seen in the case of the triplet ground state as well as for the singlet ground state, where the "zero" 
peak is however much less pronounced and the other two peaks are considerably enhanced. 

For a large level spacing (large J), only the peak at Ep is of interest. The distant side peaks would not influence 
the conductance in any way. Our calculations show that the central Kondo peak is formed when the triplet state is 
below Ep (Fig4.b), and it is not formed when only the singlet state is below Ep (Fig4.c). Such a result was expected 
because the singlet state is not degenerate, and a Kondo peak at the Fermi energy can only appear when the system 
has a degenerate occupied level (it is not exclusively a spin effect, but, more generally, it originates in the Pauli 
exclusion principle). 

In Appendix B we show that the existence of the central Kondo peak can be related with the value of the spin-spin 
correlator. The peak is absent for a strong antiferromagnetic coupling of the spins inside the dot. Such a case happens 
when the singlet state is below the Fermi energy and the triplet above. Therefore, the arguments based on the value 
of the spin-spin correlator (given in Appendix B) lead to the same conclusion as the arguments based on degeneracy 
considerations (presented in the previous paragraph). 

In this section we have focused exclusively on the description of the Kondo peaks in the DOS. However, the 
differential conductance accurately follows the DOS. In [i^, the authors prove the decoupling of Meir [s^l to give 
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FIG. 4; Plots of the local density of states (DOS) for the Kondo regime; a) eo = —0.15 and J = 0.02; b) eo = —0.11 and 
J = —0.2; c) eo — —0.05 and J — 0.2. The corresponding value of the spin-spin correlator is written on each graph. The 
vertical lines indicate the positions of the singlet and triplet levels. For all the three plots Ef = 0, F = 0.025 and T — 10~^ . 



physically correct results (qualitatively) in the Kondo regime, but the decoupling cannot be used in the mixed-valance 
regime - when charge fluctuations are mixed with spin fluctuations. This is why, in this section, we restricted the 
analysis to the spin fluctuations and the Kondo regime, when the occupancy on the dot is w 1. For the mixed- valance 
regime, the approximation made by Lacroix [i^l for the SIAM should be used instead. A generalization of [l^l was 
only done for systems with a complex noninteracting network around the impurity [46l . 0, l49l . [50l . Isil . Is^ . In the 
present paper, we added an interaction term, considerably increasing the complexity, and believe that the scheme (47j 
may no longer be applicable. The more simple decoupling used in section IV allows the analysis of charge fluctuations, 
but looses the spin fluctuations - so the two decouplings presented in this paper are complementary. 



VI. CONCLUSIONS 



The paper addresses the problem of electronic transport through a small quantum dot with a magnetic impurity. 
Our Hamiltonian model is the single impurity Anderson extended by a term of an exchange interaction between the 
electrons at QD and the magnetic impurity. The calculations have been performed by means of the Green functions 
method within the EOM scheme and presented in detail, with two different decoupling proposed, in high and low 
temperatures, respectively. 

We have studied the charge fluctuations for high temperatures, the aspect not addressed before for the case of 
QDs with magnetic impurities. The singlet and triplet conductance peaks are shown to depend on the ratio J/F (the 
exchange interaction versus the coupling strength with the leads). In particular, the transmission through the excited 
states depends strongly on the coupling to the leads, while the transmission through the ground states show just a 
broadening effect. 
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The spin-spin coupling (s • S) between the magnetic impurity and the electrons in the dot is also analyzed. The 
electron that enters the dot forms with the magnetic impurity a correlated singlet or a triplet ground state, respectively 
to the sign of the exchange parameter. The absolute value of the correlator |(s • S) \ reaches then the maximal value. 
When the excited states are shifted below the Fermi energy, the spin-spin coupling decreases. In the regime when 
the levels are dip below Ep, the spins behave like free spins and \{s- S)\ reaches its minimal value. This is due to 
an internal charge redistribution between the ground and the excited states - due to charge fluctuations: from the 
ground state to the electrode and to the excited state at QD. 

The paper also brings some new insights in the phase shift problem for the transmitted electrons. We find that 
the electronic phase shows a specific dip at the excited state position, in conditions of constant occupancy on QD. 
Similar features in the phase evolution were reported in a recent quantum interference experiment performed by 
Avinun-Kalish et al. [28l |. This experimental findings inspired Rontani [29l | to prove a generalized Fridel sum rule. 
Our calculations are an explicit demonstration that such dips in the phase evolution are due to electronic correlations 
in the system, and a redistribution of charge. The phase measurement seems to be be a very good detection method 
for such states, in a quantum interferometer device. 

In the Kondo regime, there is a competition between the spin-spin coupling with the magnetic impurity and the 
Kondo correlations with the electrons from the leads. This competition was studied before in a different system, 
namely a large isolated quantum dot, where it was found that for the exchange parameter higher than a limit value, 
the Kondo correlations are destroyed Q. For the case of a double-dot system, one can tune between a local singlet 
or a double-Kondo situation by varying different parameters of the system, like the interdot coupling [lol . [ll[. 
Our model is different, but the conclusion is in agreement: we found a limit value of the impurity-electron spin-spin 
correlator, that inhibits the formation of the Kondo resonance. For the case when both singlet and triplet levels are 
below the Fermi energy, the spins are weakly correlated and we predict three Kondo peaks in the density of states: 
one at Ep and the other two at Ep ± J, corresponding to the excitation energy. 
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APPENDIX A: SPECTRUM OF THE ISOLATED DOT 



The isolated dot is described by the Hamiltonian (1), if we put t = 0. The energy spectrum is composed from the 
empty state, the states with one and two electrons on the level eg. For single occupancy we have 







H\ n) = 


(eo + 


\j)\n) 


[Triplet] , 










(eo + 


\j)\m 


{Triplet) , 






1 U) 


+ 1 = 


(eo + 




n) + \m) 


[Triplet) , 




1 T^> 


- 1 = 






n) 


[Singlet) , 



where | ttl") is the state with the electron with the spin a =t at the state eo and the spin up at the impurity, etc. For 
double occupancy, there is a two-fold degenerate state 

H\ Ti^) = (2eo + C/)| 

F| n^) = (2eo + C/)| tm- (A2) 

For single occupancy, we see that there are two possible states. The lowest in energy will be the ground state 
of the system (singlet for positive J and triplet for negative J), and the other will be excited state. The peaks 
in transmittance correspond to the poles of the single-particle Green function, which is found with the condition 
Lj — E[i'',n ± 1) -|- E[v,n) = at u; = Ep, where E[v,n) is an energy state v from the spectrum for QD with n 
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electrons (see e.g. [53j). In our case, for Ep = 0, this is equivalent with imposing the condition E{i', n) = E{v' , n± 1), 
meaning that an energy level v from the spectrum with n electrons becomes equal with an energy level v ' for n ± 1 
electrons. This can be achieved by varying the diagonal energy eg with an applied gate potential. For the adding of 
the first electron we have the resonances at eo = 3J/4 and eo = ~J for the adding of the second electron we will 
have two resonances at eo = — C/ + J/4 and eo = —U — 3J/4. 



APPENDIX B: DISCUSSION ON THE FORMATION OF THE KONDO RESONANCE 

In order to get a deeper insight on the formation of the Kondo resonances, one can assume that, in the vicinity the 
singularity point lo ~ —a, the function Fa (Eq.l5) can be a dominating quantity in the Green function formula, and 
an expansion in 1/J^q can be done. To be more clear, we illustrate this idea first for the case J = 0, when we have the 
single-impurity Anderson Hamiltonian (within the approximations of Meir [s^ , which is similar to that of Lacroix's 
for high temperatures [i^l). We have then 



1 _.2 1 



Fo + ... . (B2) 



If the state eo is deep below the Fermi level Ep. so (no^) 0.5, it results that the coefhcient in front of Fq in the 
expression for 1/((co||cq-^)) is positive. It ensures the existence of the Kondo peak. If we look at the denominator of 
the Green function, we see that this is the condition for the — eo" line to intersect the logarithmical enhanced peak 
of F^ and therefore, to give a peak in the density of states near Ep (see Fig. I and Fig. 2, and their interpretations, 
in [47|). Of course, in the simple case J = 0, this is seen right away, the reason why we expressed the expansion of 
l/((co||cJ-|-)) is to compare the result with the same expansion for J 7^ 0, when the analytical formula of ((cot|cJ|)) is 
much more complex and not so transparent (but the expansion that we mentioned gives a result easy to discuss). 
In the case of finite J we find 

^ , ' + OM... . (B3) 



((cotI^)) l-(no,.> + (.--5) ^^^0 

The sign of the coefficient in front of Fq depends on the sign of the denominator 1 — (7104) + (s • S), because the 
numerator is always positive. In the Kondo regime, we have (?ioj.) ~ 0.5 and the denominator changes its sign for 
(s • 5"), meaning that it is positive for (s • S) > —0.5 and negative for (s • S) < —0.5. 

Three cases are presented in Fig. 4: the case a) with (s • S) ~ 0, b) (s • S) ~ 1/4, and the case c) with (s • S) ~ —3/4. 
For the case a) the coefficient in front of Fq in the expansion (17) is positive, and DOS shows the peak at Ep. In the 
case b), when the triplet is below Ep and singlet above Ep, the spin-spin average (s • S) is positive (with a maximum 
value 1/4), as the triplet state favors ferromagnetic coupling. For this case again the coefficient in front of Fq is 
positive. A different situation is encountered in the case c), when singlet is below Ep and triplet is above Ep. The 
spin-spin average is negative and, for large J, it can reach the minimal value -3/4. In this case the coefficient of Fq is 
negative, preventing the formation of the Kondo resonance. According to our calculations, the Kondo effect develops 
for {s-S) > -0.5. 
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